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Taste symmetry violations in staggered fermion formulations correlate strongly with the cut-off 
(lattice spacing) dependence in thermodynamic quantities. Better taste symmetry on the lattice 
can be achieved either by decreasing the lattice spacing and going to larger temporal extent in 
finite-temperature calculations, or by further improving the action. The highly improved stag- 
gered quark (HISQ) action offers a further degree of improvement and substantially reduces taste 
violations. We report on our studies of the 2+1 flavor QCD thermodynamics with the HISQ/tree 
action. By systematically comparing HISQ/tree, asqtad, p4 and stout calculations we quantify 
how the cut-off effects manifest themselves in different thermodynamic quantities, including the 
renormalized Polyakov loop, chiral condensate, various fluctuations and correlations of conserved 
charges. The implications for the equation of state and a comparison to the hadron resonance gas 
model are also discussed. The chiral aspects of the finite-temperature transition are discussed in 
the companion HotQCD contribution [pj]. 
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1. Introduction 

Improved staggered fermion formulations are widely used to study QCD at non-zero temper- 
ature and density [||]. Staggered fermions preserve a part of the chiral symmetry of continuum 
QCD and are inexpensive to simulate numerically After staggering the fermion field the number 
of doublers is reduced from sixteen (in four dimensions) to four. These four species of fermions, 
that we call "tastes", are inequivalent at non-zero lattice spacing, because the taste symmetry is 
broken. It is restored in the continuum limit. To reduce the number of fermions further the rooting 
procedure is used. Its validity is discussed in Ref. [ph. 

Violations of the taste symmetry at non-zero lattice spacing contribute to the cut-off depen- 
dence of different physical observables. To reduce the taste violations, smeared links, i.e. weighted 
averages of different paths on the lattice that connect neighboring points, are used in the staggered 
Dirac operator. Various improved staggered fermion formulations, like p4, asqtad, stout and HISQ 
differ in the choice of smeared gauge links. The highly improved staggered quark (HISQ) action 
was developed by the HPQCD/UKQCD collaboration and has been shown to offer a very high 
degree of improvement in the taste symmetry [f|,[5|]. This action together with the one-loop tadpole 
improved gauge action is used by the MILC collaboration ^ to generate a library of ensembles 
with several lattice spacings and light quark masses. Such a combination of the fermion and gauge 
actions was originally referred to as HISQ. In our thermodynamic studies we use a tree-level im- 
proved gauge action. Thus, to distinguish our variant of the action from the one that MILC and 
HPQCD originally used we call it the HISQ/tree action. 

In this contribution we extend our studies ^ of 2+1 flavor QCD thermodynamics with the 
HISQ/tree action, compare with the results obtained with p4, asqtad and stout and discuss the 
discretization effects in different thermodynamic quantities. 

2. Taste symmetry of different actions 

We have generated a library of zero- and finite-temperature HISQ/tree ensembles along two 
lines of constant physics (LCP): m/ = m s /5 and mi = m s /20, where mi (m s ) is the light (strange) 
quark mass. The parameters of these ensembles have been reported in [g]. For m s /5 our lattice 
spacings range from a = 0.15 fm to 0.23 fm, and for m s /20 from 0.08 fm to 0.23 fm. The statistics 
has been updated and the current status is shown in Table ||. 

Violations of the taste symmetry lead to splitting of the staggered pion spectrum into sixteen 
pion states grouped into eight multiplets, that we label in the taste basis as 75, 7075, 775, 77,-, 770, 
Yi, 7o and 1. 75 corresponds to the Goldstone pion, the lowest state whose mass approaches zero in 
the chiral limit, while the masses of the other tastes always stay finite at non-zero lattice spacing. 
In the continuum the pion spectrum becomes degenerate, and the difference in mass between tastes 
vanishes. We denote the quadratic splittings as A a = m 2 a — mL, where a is one of the above tastes 
except 75. 

In earlier work |^] we measured the pion splittings on the ensembles with the light quark mass 
of m s /5. Zero-temperature ensembles along this LCP were limited to lattices down to a = 0.15 
fm. To study the behaviour of the spectrum at smaller lattice spacings we measured A a using the 
valence light quark mass of m s /5 on zero-temperature ensembles with the light sea quark mass of 
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Table 1: The parameters of the numerical simulations: gauge coupling, strange quark mass and the number 
of time units (TU) for each run. The lattices are separated by 5 TU at zero and 10 TU at finite temperature. 
The volumes for zero-temperature runs are given in column 4. At finite temperature we used 32 3 x 8 and 
24 3 x 6 volumes. In the second column the values with error bars represent ro/a determination from the 
static quark anti-quark potential and the values marked by asterisk (*) are determined from an Allton-style 
fit. The lattice spacing is determined by using ro = 0.469 fm. 

m s /20. This is justified because quadratic splittings do not depend on the quark mass at the lowest 
order in the chiral perturbation theory, and also, as the measurements with valence HISQ on asqtad 
sea have shown, A a are sensitive to the "roughness" of configurations and not to the details of 
the action used to generate them. The results are presented in Fig. [Ij The vertical line separates 
our previously published |§] data (to the right) and the recently added calculations (to the left). 
The stout data for Ay t y s and A ri7j are shown as open symbols. We have also extracted approximate 
values of the stout Ay and Ai from the right panel of Fig. 2 in Ref. [^], and assumed for stout 
A %r5 = A^ Ts , Ayy = A^ and Ay = Ay. When plotted against a 2 , we expect the splittings to show 
some curvature for all actions, because in the scaling regime we expect them to scale with cc 2 a 2 
for asqtad and HISQ/tree, and with CC s a for stout. We found that A a for different actions can be 
fit to a polynomial form f(a 2 ) = B20: 2 + B^a A + B(,a 6 up to the lattice spacing a = 0.22 fm. After 
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performing such fits, we obtained the splittings A a as functions of a 2 for HISQ/tree, asqtad and 
stout. 

Having data for the splittings we set the Goldstone pion mass to m 75 = 140 MeV and calculate 
the root-mean-squared (RMS) pion mass as a function of the lattice spacing: 

m™ s = J ^ [m\ + m\ yi + 3m\ j5 + 3m\ y . + 3m\ 7Q + 3m\ + m\ + mfj . (2. 1) 

The results are presented in Fig. |2[ Curves show the RMS pion mass obtained from fitted values 
of A a . The thickness of each band represents the systematic error introduced by varying the end of 
the fitting interval from 0.17 to 0.22 fm. The symbols correspond to the RMS pion calculated from 
the measured pion splittings. For our estimates below we took the midpoints in each band. 

Consider lattice spacing a = 0.15 fm. On an N t = 8 lattice it corresponds to temperature 
T = 164 MeV, well in the transition region. At this a the RMS pion mass is 306 MeV for HISQ/tree, 
394 MeV for stout, and 496 MeV for asqtad. Having the same mass as for HISQ/tree requires 
a = 0. 102 fm for asqtad and a = 0.109 for stout. In other words, a HISQ/tree simulation on N t = 8 
at T = 164 MeV is comparable to an asqtad simulation on N t = 197.3/164/0.102 ~ 11.8, or a 
stout simulation on N? = 197.3/164/0.109 ~ 11.0 at the same T. Thus, we expect HISQ/tree 
Nt = 8 results to be close to stout N t = 10 and asqtad N t = 12 results, as far as the taste symmetry 
is concerned. For comparison, if one desires to have 200 MeV RMS pion at T = 164 MeV, this 
requires a = 0.089 fm for HISQ/tree, a = 0.067 fm for asqtad, and a = 0.062 fm for stout. This 
translates into the temporal extent of N t ~ 13.5, 18.0, 19.4, respectively, for these actions. The 
lattice spacings discussed above are represented with vertical lines in Fig. ||[ These estimates are 
rather crude, but are completely in line with the conclusion of ^ that a HISQ (or HISQ/tree) 
ensemble with spacing a is comparable to an asqtad ensemble with 2/3a. 

Lattice artifacts, related to taste symmetry breaking, affect masses of hadron states, and, in 
general, distort the hadron spectrum. In our simulations the masses of the pseudoscalar mesons 
m n and were used as input to constrain the LCR However, other states, e.g. vector mesons 
or baryons, are predictions and can show how well the spectrum can be reproduced at a given 
lattice spacing. In Fig. [| we present masses of p, K* , <p mesons, nucleon and £2-baryon along with 
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Figure 3: Preliminary results for masses of vector mesons (left), nucleon and Q. (middle) and the decay 
constants of pseudoscalar mesons, including the fictitious Tj sf meson (right) for HISQ/tree, m/ = m s /2Q. 
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Figure 4: The Polyakov loop (left), chiral condensate with ro (middle) and with fx (right) scale. 

a 1 extrapolations to the continuum. We have also calculated the pseudoscalar decay constants. 
Preliminary results for five ensembles are shown in the rightmost panel of Fig. pi. 



3. Taste symmetry and thermodynamic quantities 

Now we turn to the discussion of different thermodynamic quantities. In the left panel of 
Fig. |4] the renormalized Polyakov loop is shown for several actions. This quantity is associated 
with deconfinement and is not very sensitive to the chiral properties of the action. Therefore we 
observe mild effects from going to finer lattices or to actions having less taste breaking. The chiral 
condensate, shown in the middle panel of Fig. ^ on the contrary, exhibits significant differences 
for p4, asqtad and HISQ/tree on lattices with the same N t = 8. We attribute these to the amount 
of taste breaking that these actions have. P4 and asqtad, having higher degree of taste breaking, 
operate at a higher RMS pion mass than HISQ/tree and stout; therefore the transition is shifted 
to higher temperature. As shown in Sec. |2| at T = 164 MeV the asqtad RMS pion mass is about 
496 MeV, and the HISQ/tree mass is about 306 MeV. This improvement in the taste symmetry 
shifts the transition region by about 15 MeV. (Compare HISQ/tree and asqtad N t = 8 curves in 
Fig. ^[) The stout data in the middle panel of Fig. |] is plotted by converting the fg scale to ro 
scale (using fx and ro values calculated by the Budapest- Wuppertal collaboration). This allows us 
to compare the data for different actions at finite lattice spacing, when the latter is set by the same 
procedure. 

To directly compare to the results of the Budapest- Wuppertal group ^ we have also plotted 
the chiral condensate with the temperature scale set by f K (for HISQ/tree, but not asqtad) rather 
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Figure 5: The strange quark number susceptibility. 
The solid line represents the HRG model. 
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Figure 6: The light quark number susceptibility. 
The solid line represents the HRG model. 



than ro. This is shown in the right panel of Fig. As one can see, HISQ/tree N t = 6 and 8 
results collapse into one curve which is almost identical with the solid curve that represents the 
HISQ/tree continuum result extrapolated from N t = 6 and 8 using the ro scale. We attribute the 
difference between the HISQ/tree and stout continuum extrapolations to the light quark mass dif- 
ference: HISQ/tree (and asqtad) LCP uses m/ = m s /20, while stout LCP is tuned to m; = m s /27. 

We also study the strange and light quark number susceptibilities, Fig. || and [6[ These quanti- 
ties are associated with deconfmement and signal the liberation of degrees of freedom with certain 
quantum numbers. Because these quantities are related to hadron properties in the low temper- 
ature region, they are sensitive to the amount of taste breaking at finite lattice spacing. In line 
with previous observations the asqtad data at finite lattice spacing is shifted to higher temperatures, 
compared to HISQ/tree. The strangeness-baryon number correlations are shown in Fig. [7[ The 
agreement between the HRG model and the lattice is considerably improved with the HISQ/tree 
action. We conclude from Fig. [7| that this quantity is dominated by the states heavy enough, so 
that their masses are reproduced well with the HISQ/tree action even at finite lattice spacing, cor- 
responding to N t = 6 and higher. 

Finally, we discuss the effects of the taste symmetry breaking on the trace anomaly, Fig. ||. 
In the low-temperature region (T < 160 MeV) our HISQ/tree results agree well with the contin- 
uum estimate for the stout data, obtained by Budapest- Wuppertal Collaboration by averaging over 
N T = 6 and 8 or N T = 8 and 10, respectively (see [||] for details), but not with the p4 and asqtad 
results. This is expected, again, due to higher degree of taste breaking in p4 and asqtad. At high 
temperatures (T > 250 MeV) p4, asqtad and HISQ/tree results agree between themselves, but sig- 
nificantly disagree with the stout data. We note that p4, asqtad, and HISQ/tree actions improve 
the quark dispersion relation whereas the stout action (with conventional staggered fermions) does 
not. The solid line in Fig. || is a parametrization of the equation of state derived in [^] based on the 
hadron resonance gas model and lattice data. 

4. Conclusion 



Taste symmetry breaking at finite lattice spacing is one of the major cut-off effects in staggered 
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Figure 7: The S-B correlations. Figure 8: The trace anomaly. See details in the text. 



actions. By comparing several staggered actions with a different degree of improvement and at 
different lattice spacings, we studied the cut-off dependence of several thermodynamic quantities. 
As expected, the quantites that are either directly related to the chiral symmetry (such as the chiral 
condensate) or are sensitive to the hadron spectrum (for instance, quark number susceptibilities and 
the interaction measure) are affected by taste symmetry breaking of the staggered fermion actions. 
The HISQ action, having the highest degree of improvement, exhibits the smallest cut-off effects, 
compared to p4, asqtad and stout. 
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